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We investigate bosonic sectors of supersymmetric field theories. We consider superpotentials 
described by one and by two real scalar fields, and we show how the equations of motion can be 
factorized into a family of first order Bogomol'nyi equations, so that all the topological defects are 
of the Bogomornyi-Prasad-Sommerfield type. We examine explicit models, that engender the Zjv 
symmetry, and we identify all the topological sectors, illustrating their integrability. 
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PACS numbers: 11.27.+d, ll.30.Pb, ll.30.Er 

Domain walls are defect structures that appear in diverse branches of physics. They usually live in three spatial 
dimensions as bidimensional objects that arise in systems described by potentials that contain at least two isolated 
degenerate minima. They involve energy scales as different as the ones that appear in Condensed Matter Q and in 
Cosmology ||]. 

A lot of attention has been drawn recently to domain walls in field theories, in models that have been investigated 
under several distinct motivations |]3|-|l5|| . A very specific motivation concerns the presence of domain walls arising 
in between non zero vacuum expectation values of scalar fields in supergravity ||^,^. Another line deals with the 
formation of defects inside domain walls ■ A great deal of attention has also been drawn to SU (N) supersymmetric 
gluodynamics, where nonperturbative effects give rise to gluino condensates that may form according to a set of N 
isolated degenerate chirally asymmetric vacua, from where domain walls spring interpolating between pairs of vacua 

0- . . . . . 

The interest in domain walls in general widens because of the interplay between Field Theory and the low energy 

world volume dynamics of branes in String Theory |l^~[ig|]. In the case of intersection of defects, in particular in 
Ref. Jl2| some aspects of wall junctions have been investigated when the discrete symmetry is the Zn symmetry, in 
systems described by a single complex field, with the superpotential 
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The second work in Ref. 
be cast to the form 



p2(| has shown that the tensions of the topological defects that appear in these systems can 
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where [N/2] is the biggest integer not bigger than N/2 itself. In Ref. jsj the same result was obtained, in an 
investigation concerned with exact integrability, due to the presence of infinitely many conserved currents - see also 
Ref. [gOj for other investigations, in particular on the explicit form of the BPS solutions at large N. 

In the present work we examine the bosonic portion of supersymmetric theories, similar to the above models, that 
are of general interest to supersymmetry. Our investigations follow the lines of the former work |p^ , and bring new 
results on the presence of domain defects in models described by one and two real scalar fields. We start investigating 
systems with a single real scalar field. We write the Lagrangian density in the standard form 
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W = W{(j)) is the superpotential, and = dW/dcf). We search for defect structures, for static solutions of the 
equation of motion. We suppose the static solutions depend only on x, on a single spatial coordinate. The equation of 
motion becomes cP'(t>/dx^ = W^W^^. We examine the energy of the static configurations, and we associate with the 
second order ordinary differential equation of motion the two first order ordinary differential equations d(j)/dx — ±W^. 
These first order equations are Bogomol'nyi equations and their solutions are named BPS states, which saturate the 
lower bound in energy. 
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We see that solutions of the first order equations solve the second order equation of motion, and the proof follows by 
direct differentiation of the first order equations. We also see that one associates with the second order equation two 
first order equations. This fact suggests the possibility of factorizing the equation of motion into first order Bogomol'nyi 
equations. This is indeed the case, and we demonstrate this property by examining the ratio R{4>) = {d(f> / dx) / {dW / d(j)) . 
We use this expression to write 
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We notice that dW'^/dx = 2{d(j) / dx)W^W^^ and that 

where we have used the equation of motion. These results show that the quantity S{(f)) = — (dip/dx)"^ does not 
depend on x when (j){x) solves the equation of motion. For those 0(a;) we have that limx^-oo 4'{x) ~^ : with 
being a vacuum state, and \im.x^-oa{d(j)/ dx) = 0. Also, lima:^_oo = 0, since the vacuum states are extrema of 
the superpotential. Thus we get that S{(j)) vanishes, and this allows writing R{cj)) = ±1, which gives the first order 
Bogomol'nyi equations. This result shows that the second order equation of motion is completely equivalent to the 
two first order Bogomol'nyi equations. A direct consequence of this result is that in such models all the topological 
solutions are of the BPS type, that is, these models do not admit the presence of non-BPS states. The first order 
Bogomol'nyi equations can be readily integrated, so our result offers a new way of showing exact integrability. 

We illustrate our result with some examples. There are several models of systems of a single real scalar field. We 
work with dimensionless quantities, and we consider the superpotential 

7^ = 0,1,... (6) 

n + 1 n + 6 

Widely known examples are described by n = 0, 1, and identify the (f)^ and the 0^ models, respectively. These models 
engender the Z2 symmetry, and for n = there are two asymmetric vacua, = —1 and = 1. For n = 1,2, ... 
there are three vacua, the former two and another one, symmetric, = 0, that appear with multiplicity n. For n = 
the first order equations give the standard defects, 4){x) = ±tanh(x). These BPS solutions connect the two distinct 
asymmetric vacua. For n — 1 the solutions are such that (jy^ix) = (1/2)[1 ± tanh(a;)]. These BPS states connect the 
symmetric vacuum wq = to the asymmetric ones, w^'^ = ±1. The models defined by the above superpotential are 
exactly solvable; the solutions can be known implicitly for n = 2,3, ... The tensions are: for n — 0, t'^ — 4/3, and for 
n = 1,2, = 2/{n + l){n + 3). They show that the width of the defect increases for increasing n. 

The main difficulty to generalize the former result is that it is unclear how to write something like the ratio R{4>) 
if we deal with two or more fields. However, we have found an interesting possibility, that appears in the case of two 
real scalar fields, when the superpotential is harmonic on the two fields. The issue here is that for Wicf), x) harmonic 
we can work with a complex superpotential, written as W{(p), in terms of the complex field (p — (f> + ix ^ see Ref. [ p^ . 
In this case the ratio R{(p) generalizes to R{f), in the form 

R{^) = — ? (7) 
W'{ip) dx 

The superpotential W{lp) is holomorphic, and we consider systems that are described by the Lagrangian density 

C^\d^-^d'^^-\\W'{^)\^ (8) 

where W'{if) — dW/dip. The specific form of the potential shows that the critical values of the superpotential are 
minima of the potential. We search for defect structures, for static solutions of the equation of motion. We suppose 
the static solutions depend only on x, so the equation of motion becomes 

l^ = W'[^)W^) (9) 

We examine the energy of the static configurations to show that the second order ordinary differential equation is 
associated with a family of first order ordinary differential equations. They are 
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where ^ is some real constant. These first order equations are Bogomol'nyi equations. As usual, solutions of these 
first order equations also solve the second order equation of motion. The proof goes as follows: we differentiate the 
first order Eq. to get 



dx-^ 



dx 



W'{ip)e 



■ dJp 



= W"(ip) ^ e''^ = W"(ip) W'(ip) 
dx 



(11) 



In the above equation, in the last equality we have used the first order Bogomol'nyi Eq. (|T^) once again. 

The above models are somehow similar to the former models, described by a single real scalar field. Thus they may 
also lead to factorization of the equations of motion into first order Bogomol'nyi equations. We make this reasoning 
unambiguous by concentrating on the ratio R{(p) defined by Eq. (^. We can write 
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where we have used the equation of motion. We use the above Eqs. (^3|) and (|TJ) to see that now the quantity 
S{p) = \ W'{p)\'^ — \dp/dx\'^ is such that dS{p)/dx = 0, that is, it does not depend on x when p{x) solves the equation 
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with V a vacuum state, and lim^;- 
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of motion. For those p{x) we have that lima,- 
Also, lim2;_>_oo W {p) = 0, because the vacuum states are critical points of the superpotential. Thus we get that S{p) 
vanishes, so that R{p) does not depend on x. Also, \dp/dx\ = \W' {p)\ and this allows writing \R{tp)\ — 1. Thus we 
get that R{ip) = e~'^ for some constant ^, which leads to the first order Bogomol'nyi equations. 

This result shows that the second order equation of motion is equivalent to the first order Bogomol'nyi equations. 
This is valid with the boundary conditions lima;^_oo ^{x) — > , Ximx^-ooidf/dx) = 0, as required by the topological 
solutions. A direct consequence of this result is that in such models all the topological solutions are of the BPS type, 
that is, these models do not support non-BPS states. Furthermore, we recall that in systems of real scalar fields, one 
may occasionally run into solutions that engender no topological feature. They are named nontopological solutions, 
and have been found for instance in Ref. We first notice that in systems of real scalar fields described by 

some superpotential, the nontopological solutions cannot appear as solutions of the first order Bogomol'nyi equations, 
because they should have zero energy, and this is the energy of the vacuum states. We add this to our former result 
to obtain another result, that in systems of two real scalar fields, when the superpotential is harmonic there is no 
room for non topological solutions. 

To illustrate these results let us consider the superpotential 



W^{p) 



1 



1 



1 



N - 



1 



It is defined by -/V = 2, 3, ... and hy n — Q, 1, ... The potential is given by 

v^\{^pT{i~p^) (i^p^) 



(15) 



(16) 



The models engender the Zm symmetry, irrespective of the specific value of n, and the limit n leads us back to 
the model described by Eq. (ph. Below we obtain the tension associated with every BPS solution in a closed form. 
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The critical values of WJ^{(p) can be readily obtained. For n = they are v'^ = exp(i^^), where = 27r(fc/iV) 
and k = 1,2, A*". For n 7^ 0, we have to add to this set of extrema the special value vq = 0, that appears with 
multiplicity n. Thus, while in the case n = we get N asymmetric phases, degenerate, in the case n 7^ we have 
to include another phase, that breaks no symmetry of the original system. This includes the possibility of a chirally 
symmetric phase in SU{N) supersymmetric gluodynamics. 

We investigate the tension of the defect structures. Since the defects solve first order Bogomol'nyi equations, their 
tensions are obtained as: for n ^ 0, for solutions that connect the origin vo to any of the minima we can write 
t%Q = \ WJ^{1) I, that define the radial sectors. The result is, for N arbitrary, and for n = 1,2, ... 



N 

(n + l)(iV + n+l) 



We set = 2 to get — 2/(n + l)(n + 3), which reproduces former result, obtained below Eq. (y). 

The other sectors represent defects that connect vacua in the unit circle. Here we have to calculate the quantity 
\W^{v^) — W^{v'^)\. The result can be cast to the form 



(n + 1) — 
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where k = 1,2, [N/2] classifies the minima on the circle that the defect connects: fc = 1 for first neighbours, k — 2 
for second neighbours, and so forth. We set n = to get back to the former result, Eq. (||). 

We use the superpotential (|2^) to write the family of Bogomol'nyi equations in the form, in terms of the real fields 
(j) and X, 



dx 
dx 



P:;^cosC + Q^sin^ (19) 
Q^cos^-P^sin^ (20) 



Here we have introduced the quantities 



and 



= P„(l-PAr) + Q„Qw (21) 



= Qnil - Pn) - PuQn (22) 

which are defined in terms of the functions 

Pfc = 7^e (0 - ^X)^ Qk^Im{dp-ixf (23) 

We illustrate the general situation considering particular values of N and n. The case n = is special, because it 
does not admit the origin in configuration space as a vacuum state. We consider n = and N even, N = 2j, j = 1,2, ... . 
We examine the Bogomol'nyi equations for solutions that connect the vacua (±1,0), which are farthest neighbours 
in the set of iV = 2j minima. In this case the equations demand that x = 0, so there are only one-field solutions, 
obeying d(f)/dx = ±(1 — (f)^^). In the case j = 1 we get to the standard tanh defects that we have already described. 
For j = 2, 3, ... we have other solutions, that can be known implicitly. They behave as (t>{x) = x — x'^^'^^ /{2j + 1) in 
the limit a; — > 0. They are similar to the standard defect that we have already found. However, they are thinner and 
then more energetic as j increases, as it can be readily verified. We use the tension result to get t^j j = 4j/(2j + 1). 
This tension increases from 4/3 for j — 1 up to the finite value 2 in the limit j 00. There are other solutions, that 
connect other vacua in the unit circle. 

In particular we consider the case N = 2,n = 0. The superpotential is — (p^/i + (px^ . It is instructive to see 
this model with another superpotential, more general, W{(t), x) = — 0'^/3 — ripx^- The model with = 2, n = is 
obtained for r = — 1, the value that makes the superpotential harmonic. For r G (0, 1/2), however, the more general 
model maps the anysotropic XY model 0]: there are one-field solutions describing Ising walls, and two-field solutions 
describing Bloch walls - see the last work in Ref. ||l2|. For r — —1, for non vanishing x one must have (jP — 1 + x^/3, 
so there is no two-field solution connecting the minima (±1,0) for r = —1. The lesson we learn here is that the 
condition W^^ + = certainly simplifies the model, favoring its integrability. 

Similar behavior is also present in the case N = 2, for n arbitrary. Here all the vacuum states are in the axis. 
The first order Bogomol'nyi equations only support topological solutions for x = 0, irrespective of the value of n. The 
Bogomol'nyi equations effectively reduce to the form 
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which supports no nontrivial two-field solutions. They are the first order equations that appear in the superpotential 
of Eq. (^). The result is that for N = 2 and n arbitrary the systems reduce to systems of just one field. The tensions 
of the defects for TV = 2 are: for n = 0, ^ = 4/3, and for n ^ 0, t^^ = 2/{n + l)(n + 3), as noted below Eq.(|7|). 

We examine ^ to see that it vanishes for some values of N, n, and k. The condition n = — 1 is special, because 
now the models only support N radial sectors, that describe solutions connecting (0, 0) to any of the N vacua in 
the unit circle. For n = N — 1 the BPS sectors effectively reduce to the sector described by just one field, and the 
Bogomol'nyi equations become 

= ±0^-1(1-0^) (25) 

ax 

The corresponding tensions are i^g^ ~ 1/2N. The case N = 2 reproduces the case n = 1 of the superpotential in 
Eq. (^. Other values of N and n give rise to radial sectors, and to sectors that connect vacua in the unit circle. An 
example of this is given by iV = 6 and n = 2. This model supports six radial sectors, with tension <| q = 2/9 and six 
other sectors, connecting first neighbours in the unit circle, with tension i ^ 4/9, twice the tension of the radial 
solutions. 

Evidently, the new idea of factorizing equations of motion into first order Bogomol'nyi equations is of direct interest 
to supersymmetry. It can be used in diverse applications, in particular to investigate the presence of BPS walls 
interpolating between distinct pairs of vacua, and exact integrability. The complete factorization of equations of 
motion in systems described by W{<j)^ x) harmonic, or by W(}p) holomorphic, ensures BPS feature to domain walls 
that spring in such systems. Because of the associated BPS character, these domain walls partially preserve the 
supersymmetry. 

The superpotential W[cf), x) in general depends on a set of parameters, and these parameters can be adjusted to 
make W{4>, x) harmonic or not. An example of this can be found for instance in the second paragraph below Eq. (|2^), 
in the model used to map the anisotropic XY model, which supports Ising and Bloch walls - in the case r — —1 
the superpotential there considered becomes harmonic, but the model does not support Bloch wall anymore. On the 
other hand, harmonic and non-harmonic superpotetials of the form W{(j)^x) may be used to distinguish two types of 
supersymmetry, one with M ^ 1 and the other, enhanced to the case M = 2. This behavior has been examined for 
instance in Ref. p2| . In that work, we can also find Wess-Zumino models with two complex fields, and there it was 
shown that the system supports topological sectors of the non-BPS type. Thus, the non-BPS solutions can be used 
as examples that show that our result does not hold in Wess-Zumino models with two or more complex fields. 
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